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The critical effective potential is the nonperturbative part of the effective action at a phase transition. 
It equals the scale invariant effective average potential and can be calculated from the renormaliza- 
tion group flow of the effective average action. In some cases this requires only the solution of an 
ordinary differential equation without actually simulating the renormalization group flow. Here the 
Ising model is examined beyond leading order and with full field dependent effective potential. 
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The effectivepotential in statistical physics and quantum 
field theory [l|, y includes fluctuations of all wavelenghts 
as opposed to the microscopic potential. It is therefore 
relevant for the computation of macroscopic properties. 

The usual way to access the effective potential is to start 
vifith the potential given by the microscopic theory and 
follow the renormaUzation group flow to the infrared . 
As the potential then depends on fleld variables and the 
scale, this requires solving a partial differential equation. 

However, an especially interesting case is a second order 
phase transition where the correlation length diverges. 
The system becomes therefore scale invariant and the 
renormalization group flow ends up in a fix point, the 
effective critical potential. Its calculation requires only 
the solution of an ordinary differential equation |j, |y, 



Effective average action 

The effective average action ( see e.g. and references 
therein ) Tk interpolates between the microscopical or 
classical action S ~ Poo and the effective action F = 
Fq. By construction it includes only fluctuations with 
momenta larger than k and this transition is described 
by an exact renormalization group flow equation: 



(2) 

Here FL is the two point function and Rk is an arbitrary 



momentum cutoff with the properties Rk 
Rk oo for fc — > cxD and Rk{q^) > for 



for k 
0. 



0, 



We want to consider a particular interesting and common 
model, the three dimensional Ising model. The 0(1) sym- 
metry of this theory can be exploited for a lowest order 
derivative expansion of the effective average action: 



in terms of the most general form of the nonperturbative 
effective average potential Uk{p) and a field independent 
but scale dependent wave function renormalization Zk- 

When inserted into the flow equation for the effective 
average action this ansatz gives the flow of the effective 
average potential. 



UkUk - 2 J (2^)d Ul+2pU'^+q-^Zk + Rk{q'^} 



As the potential is only determined up to a constant, it 
is convenient to consider its first derivative instead. The 
corresponding flow equation is given by straightforward 
differentiation: 

"kUk - 2 J (2¥F {Ul+2pU'^+q^Z^+R^(q^)y^ 

Switching to dimensionless quantities u' — U'j,/ {Zkk"^) 
and choosing a linear cutoff function Rk — Zk{k^ — 
q^)Q{k'^ — q^) 0, ^3 the flow equation becomes ex- 
plicitely scale invariant: 

dtu'{x) ^ (-2 + ri)u' + {d-2 + r])xu" 

_ d+2-i) 3m"+2xm'" 
^'^d.d{d+2) {l+u'+2xu"y-' 

Here t = Ln(A:), x ^ Zkk^-'^p, v^^ ^ 2'^+^Tr'^/^r{d/2) 
and rj denotes the anomalous dimension r/ = dtl^nZk- 

The latter one can either be taken self consistently from 

y — ~r'^ (1+2kA)-4 

( where K(t) denotes the running minimum ( — ) 

and A = u"{k), u^^^ = u"'{k) ) or from the quite accurate 
approximation 

which will be used in this paper. 

For comparison with literature we also use a sharp cutoff, 
where in leading order 



F, 



/ d''x(Uk{p) + Zk\^^'^{x)^M< 



dtu' = -2u' + {d~ 2)xu" - 2vd 



3u" -{-2xu"' 
l-\-u' -\-2xu" 
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Scaling solutions 

The scale invariance of some system at a second order 
phase transition manifests itself in a fix point in the 
renormalization group flow of the dimensionless, scale in- 
variant flow equations. We are thus interested in scaling 
solutions u' such that 



— dfu' — acu' - 
where 

_ d(d+2) 2-1) 



bcxu" 



, 3u"+2xu"' 
'(l+u'+2xu"Y 



4va d+2-Tj 



, _ d(d+2) -d+2-ri 
~ 4va d+2~ij 



This can be rewritten canonically in terms of first order 
ordinary differential equations for u'{x) and u"{x): 

d^u' = u" 

d^u" = 27 ((a"' + bxu"){l + u' + 2xu"f + 3ii") 

Given initial values Uq, Uq at some point xo the differ- 
ential equations yield a local solution m', u" in a vicinity 
around xq. The space of such solutions has dimension 
two, its elements being uniquely defined by Uq, Wq. 

However, we are interested only in physical solutions 
which are continuations of the local solutions to the whole 
positive real axis. In other words, we must find a point 
xq > and initial values Wq, Ug such that u' is a smooth 
function for all x > 0. 

Such a search is nontrivial, though, even numerically. As 
soon as the space of local solutions has a dimension larger 
than one, global solutions can be missed even if one is 
already close to the correct initial values. Fortunately, 
by using the constraint of requiring a global solution we 
are able to eliminate one of the two integration constants: 



If dtu' — Vx, this holds true especially for xq = 0: 

= dtu' \xa = acug + cj^iji 
And hence 

In the analogous calculation for the sharp cutoff we find: 

Q^dtu' ^ acu' + bcxu" + c ?""+^o"'", 
where a = v^^, b = f^^(l — d/2) such that 

dxu' = u" 

d^u" = ^ [{au' -f W)(l + u' + 2xu") + 3u") 



Furthermore, we have additional constraints on u'^. the 
physically interesting u' has a single zero and is positive 
for large x. Thus u'q < 0. On the other hand, the denom- 
inator 1 + u' + 2xu" must be nonzero for all x. Therefore 
u'q > — 1 . Now the global solution can be found be vary- 
ing only u'q, in the range —I < u'q < 0, until the above 
requirements are met. This can be done numerically ef- 
ficiently. 

For any initial u'q the solution u' extends to some Xmax . 
For too small u'q the flow diverges to positive infinity. On 
the other hand, for large initial values the flow crosses 
zero a second time. These two cases are separated by a 
critical initial value which leads to the physical solution 
we are interested in. The plot of x^ax a-s a function of 
u'q is then sharply peaked around this critical value: 
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FIG. I: 

Xmax{u'o) foT the linear cutoff (a), 
77 = (b) and the sharp cutoff (c) 



In this way we find: 

linear cutoff, 77 = 1/27 (a) 
linear cutoff, r] — (h) 
sharp cutoff, = (c) 

In perfect agreement with 

reference [11] ( 77 = ) 
reference [5j ( 77 = ) 



-0.16902863438. 
-0.18606424944., 
-0.46153372007., 



Ai*=-0. 1860642. 
cr=-0.46153372... 
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Effective scaling potentials 



The critical exponent 6 



We can now plot the solutions u' for the above initial 
values: 
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FIG. II: 
u'{x) for the linear cutoff (a), 
1] = (b) and the sharp cutoff 
Compare plot (b) to fig. 2 in 



Due to the numerically limited accuracy of the initial 
value and the integration we are able to track these solu- 
tions only up to X w 0.12 in case of the linear cutoff and 
a: « 0.2 in case of the sharp cutoff. 

The second quantity to be read off from the critical ef- 
fective potential is its zero k where u'{k) = 0: 



linear cutoff, ij = 1/27 (a) 
linear cutoff, = (b) 
sharp cutoff, i] = (c) 



K=0.03060106819. 
K=0.03064764922. 
K=0.04711134650. 



The value k=0. 03060... is in good agreement with the fix 
point value k*=:0. 03053... from the dynamical method. 

As one can see from the different solutions for the linear 
and the sharp cutoff, the effective potential is not cutoff 
independent. Interestingly however, the location of the 
minimum hardly depends on the anomalous dimension in 
case of the linear cutoff. Moreover, it is also only weakly 
dependent on the initial value u'q. 



That the effective scaling potential in fact describes 
critical phenomena can be seen from the power law 
u' ~ x^'^"^^/^ for large x. 

The critical exponents S and rj are related by 

c _ d+2-r, 

Hence we expect 



lim 1 + 2 



d\n(u') _ 



dln{x) 
in case 77 = and 



= 5 



lim 1 + 2 



d\n{u') 



67 
14 



dln{x) 
in case 77 = 1/27. 

These asymptotic values are plotted as horizontal lines: 
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FIG. Ill: 

1 + 2dhiu' /d\nx for the linear cutoff (a), 
rj — (b) and the sharp cutoff (c) 
approaching the critical exponent 5. 

The asymptotic behaviour is as expected in all three 
cases. However, if one wanted to calculate S from the po- 
tential its precision would suffer from the limited range 
of integration. 



Comparison to the dynamical method 

In order to assess the advantage in computing complexity 
of the method presented, we also take a look at the usual 
dynamical method 0, 0| . In this case one starts with 
an initial potential = Aa(x — ka) at some ultraviolet 
scale A with an initial expectation value of the fields ka 
and some quartic coupling constant Aa. 

With this initial potential the partial differential equa- 
tion for dtu'{x) is computed on a discretized grid. The 
renormalization group flow for t — > — cxd then either leads 
to the symmetric phase ( k — > ) or to the phase with 
spontaneous symmetry breaking [ k ^ oo ). 

These phases are separated by a second order phase tran- 
sition. Correspondingly, there is a critical ka = Hcr- Fine 
tuning ka this value can be found. The renormalization 
group flow near the phase transition will then come close 
to the fix point and the scaling potential can be read out. 

The fine tuning of ka is now replaced by fine tuning /q. 
In effect, the solution of a partial differential equation is 
replaced by the solution of an ordinary differential equa- 
tion. This is much simpler and may reduce computing 
time by several orders of magnitude. 



Conclusion and outlook 

The method allows the calculation of the effective scaling 
potential with very little effort. The integrable range 
seems to depend only on numerical precision and could 
maybe be improved by the use of high accuracy floating 
point libraries. 

The effective scaling potential could also be used as an 
initial value in a dynamical simulation, allowing the de- 
termination of critical exponents depcndend on small de- 
viations from the critical point, while still overcoming 
the time consuming necessity for a dynamical fix point 
search. 
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